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A variable speed of light (VSL) cosmology is described in which the causal mechanism of gener¬ 
ating primordial perturbations is achieved by varying the speed of light in a primordial epoch. This 
yields an alternative to inflation for explaining the formation of the cosmic microwave background 
(CMB) and the large scale structure (LSS) of the universe. We make use of the SAT formalism 
to identify signatures of primordial nonlinear fluctuations, and this allows the VSL model to be 
distinguished from inflationary models. In particular, we find that the parameter /nl = 5 in the 
variable speed of light cosmology. The value of the parameter <?nl evolves during the primordial era 
and shows a running behavior. 
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I. INTRODUCTION 


Inflationary cosmology is successful in explaining the 
CMB observational data. However, it is not the unique 
paradigm of the very early universe. For example, we 
have bouncing cosmology, string gas cosmology, and the 
variable speed of light (VSL) cosmology [l|. We refer 
to refs. [3H6| for comprehensive reviews on bouncing cos¬ 
mologies and to !§} for reviews on string gas cosmology. 
It is necessary to obtain an observational signature that 
can distinguish between the different competing models 
of the early universe. If primordial B-mode polarization 
gravitational waves can be detected, then the measured 
ratio r = tensor/scalar will help to distinguish between 
the models, for this requires a mechanism to produce a 
stretching of both scalar perturbation density and gravi¬ 
tational waves to superhorizon scales [l| , so that they can 
be detected in the CMB. However, a measurement of the 
tensor mode tilt nt is needed to give a definitive test of 
the models. A joint analysis of the BICEP2/Keck Array 
observational data found strong evidence for foregound 
dust and no statistical significant evidence for gravita¬ 
tional wave tensor modes. An upper limit was obtained 
for the tensor-to-scalar mode ratio r < 0.12 [o] 


Another important test of competing models is the 
magnitude of nonlinear perturbations of the scalar and 
tensor power spectra. In the following, we will calculate 
the magnitude of the nonlinear perturbations in a ver¬ 
sion of VSL U, using the so-called SAf formalism. We 
will find that the nonlinear perturbations are of order 
unity, and the non-gaussian parameter /nl and its run¬ 
ning behavior g^L are consistent with the Planck Mission 

data (Toj . 


II. THE VSL MODEL 


The model is formulated using a Friedmann-Lemaitre- 
Robertson-Walker FLRW metric: 


ds 2 = c 2 dt 2 - 


dr 2 


1 - Kr 2 


r 2 (d0 2 + sin 2 Odcf 2 ) 


, (1) 


with K = 0, +1, —1 for flat, closed and open models. The 
metric has the group symmetry 0(3) x R with a preferred 
proper comoving time t. 

Setting the cosmological constant A = 0, the variation 
of the action with respect to g^ yields the field equa¬ 
tion [l|: 


G 


fllS 


8t rG^ , V“V a $ 

$ 1 ^ + $ 9 ^ $ ’ 


( 2 ) 


where G^ — \g tlv R- V p is the covariant derivative 

with respect to the metric g^ u and $(x) = c 4 (x). We also 
obtain the field equation for the scalar “seed” field <f >: 1 




dV{4>) 

d(j) 


( 4 ) 


where V (</>) is a potential. In the action S, we have made 
the speed of light c(x) a dynamical degree of freedom. 
The energy-momentum tensor T Ml/ satisfies the conser¬ 
vation law: 


= 0 . 


( 5 ) 


1 From the action principle [jj, we also have the field equation: 


_ aVFQ/y) = 

&ll>v 


(3) 


where 0^ is a vector field, B^ v = du,i\) v — and W(i^^) 

is a potential. For (010^ |0) ^ 0 the Lorentz group SO( 3,1) is 
spontaneously broken down to 0(3). 
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We postulate that fractions of seconds after the big 
bang c undergoes a large superluminal increase c > eg 
where Cq is the present measured speed of light. This 
solves the horizon and flatness problems [l]. Moreover, 
the primordial quantum scalar matter fluctuations and 
gravitational wavelengths X s and A„,, respectively, are 
stretched superhorizon for c > Co and frozen in as clas¬ 
sical fluctuation imprints on the CMB: 

c c 

Xs — 5 X w — , ( 6 ) 

v s v w 


where v s and v w denote the constant values of the scalar 
and gravitational wave frequencies, respectively. 

The Friedmann equations in our model are given by 


H 2 


Kc 2 




( 7 ) 


and 


III. NON-GAUSSIANITIES IN THE VSL 
COSMOLOGY 

In this section, we first briefly review the standard pro¬ 
cess of calculating nonlinear perturbations by virtue of 
the so-called SAf formalism. Within the context of the 
VSL cosmology, primordial fluctuations are seeded by the 
scalar field (f> and, therefore, the mechanism of generating 
primordial curvature perturbation in the VSL cosmology 
is analogues to the curvaton scenario emu. Then, with 
the assumption that the end of the epoch c ^$> c o when 
c = Co is a uniform total density slice, we generalize the 
SAT formalism jl5l - [l9j to be applicable to our case. Af¬ 
terwards, we analyze the non-gaussianities of primordial 
curvature perturbations in VSL cosmology. 


A. The generalized 5Af Formalism 


a 

a 



( 8 ) 


where H = a/a and p = pm + pip + pp>- 

The scalar field <p that describes the structure growth 
“seed” held satisfies the Held equation: 


4> + 3H<j) + c 2 V{cj>) = 0. (9) 


We assume that V{<j>) =0 and obtain the solution: 

j> = VnB^j , (io) 


where B is a constant and a* is a reference value for 
a. For large p the kinetic contribution to pp, oc \{4>) 2 
will dominate the matter densities pm and pp, and the 
Friedmann equation 0 becomes 


H z 


Kc 2 



Substituting the solution (flUl) . we get 


H 2 


Kc 2 



( 11 ) 


( 12 ) 


Because the held <p dominates in the early universe, we 
can neglect the spatial curvature for c = cq and we obtain 
the approximate solution given by 

a(t) = a*(3Bt) 1/3 , H(t) = ^. (13) 

The equation of state for a massless scalar held gives 
for the exponent n in a(t) oc t n the value n = 2/3(1 + 
w), so for n — 1/3 we get w = 1. A calculation of the 
power spectra and spectral indices (tilts) for the scalar 
perturbation density fluctuations and gravitational wave 
tensor fluctuations yields: n s = 0.96 and the red tilt 
n t = —0.04 [J. If we adopt the tensor/scalar ratio bound 
r < 0.12 then we have \r/n t \ < 3. 


The SAf formalism is available upon two conditions: 
hrstly, the universe is isotropic and homogenous at ex¬ 
tremely large scales that can accommodate a large num¬ 
ber of causally isolated regions; secondly, the perturba¬ 
tion is frozen after the horizon exit. These two conditions 
are easily met in inhationary cosmology. However, for al¬ 
ternatives to inflation, if one can argue based on them, 
the SAf formalism is still available. The hrst condition 
is in fact related to the horizon problem and, hence can 
be satished in quite a number of alternative models. For 
the second condition, it is not satished in matter bounce 
cosmology and marginally satished in ekpyrotic cosmol¬ 
ogy. In VSL models, the universe is still expanding with 
an effective equation of state parameter w = 1, and as 
a result, the dominant perturbation grows as a logarith¬ 
mic function. In this case, the second condition is also 
satished and from this argument the SAf formalism is 
available as well in VSL. 

On superhorizon scales, if we assume negligible inter¬ 
action between primordial scalar held tp and other matter 
helds, the component curvature perturbation on the uni¬ 
form density slice can be identihed as 


C i( x ) 



dpi 

Pi + Pi{pi) ’ 


(14) 


through the SAf formalism, where the subscript i repre¬ 
sents the component matter held in the universe. 

Having the component curvature perturbation in mind, 
we can calculate the curvature perturbation on the uni¬ 
form density slice. In analogy with the curvaton mecha¬ 
nism, it is natural to choose this slice at the end moment 
of c Co period when c = Co- On the uniform density 
slice we have 


Pr+ P<p = P , (15) 

where p r and pp, represent the radiation and scalar held 
densities, respectively. For the primordial scalar held, 
the background dynamics behaves as a stiff fluid with an 
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effective equation of state w<f, = 1. Then, following Eq. 
m, we can obtain [ 20 | : 


Moreover, the two point correlation function is related to 
the regular power spectrum through the relation: 


n - n P 4 (<h-C) 

rr i 


p<!> = c) . 


(16) 


As a consequence, the curvature perturbation ( can be 
derived on the uniform density slice through the relation: 


C = r(d„ 

where 

_ _ 3 
2 +12^ 


(17) 

(18) 


Here, we have introduced the transfer efficiency coeffi¬ 
cient r and = p^/ptot is the density parameter for 
the primordial scalar field. Note that, in a generic case, 
f = 3(1 + w ( j,)fl ( j,/[ 4+ (3 Wc/, — 1)12^], and we have applied 
the relation w$ = 1. In the VSL cosmology, we have 
Vtff, = 1 , for the scalar field dominates over other mat¬ 
ter fields at the end of the primordial era. This implies 
that the transfer from the scalar field fluctuations into 
primordial curvature perturbations is very efficient and 
thus, r = 1 . 

Within the local ansatz of the curvature perturbation, 
one can expand its form order by order as follows, 


C(x) = CiO) + C2O) + C3O) + G(C4) 

= CiO) + it/nlCiOe) + ^9nlCi ( a: ) + 0(C4), (19) 

where £1 is the fluctuation of the Gaussian distribution, 
and £„ are the non-Gaussian fluctuations. From the 
above expansion, we obtain 


nm = 


(23) 


B. Primordial perturbations and non-Gaussianities 
in the VSL cosmology 

Following the discussion in the previous subsection, 
one can find that the dynamics of primordial curvature 
perturbation would be identified by the variation of the 
generalized VSL e-folding number by making use of the 
local ansatz. To calculate such a variation, we need to 
know the evolution of the background universe. The 
scalar field satisfies the Klein-Gordon equation: 

$ + 3H0 = 0. (24) 

In this phase, a oc 2 1 / 3 and, therefore, the Hubble param¬ 
eter is expressed as 


H ® = Yf (25) 

According to the convention of ref. |l| • the primordial 
scalar field <fi is dimensionless. Substituting Eq. (l25l) 
into the background equation of motion (l24l) yields the 
following solution: 


(j)[t) = </>/+ 7 In 



4>{t) = j, 


(26) 


/NL=g^2, 5 nl = y^3, (20) 

for the local type. To be general, one can relate the non¬ 
linearity parameters to the bispectrum and the trispec¬ 
trum via: 

-B(k 1; k 2 ,k 3 ) = ^f N L[P(ki)P(k 2 ) + 2perm] , 

54 

T(ki,k 2 ,k 3 ,k 4 ) = — g N L[P(ki)P(k 2 )P(k 3 ) + 3perm] 
25 

+ tnl [-P(fci)P(fc 2 )P(|ki + k 3 |) + llperrn], 

( 21 ) 

where these spectra are associated with the correlation 
functions as follows: 

2 

<C(k 1 )C(k 2 )) = (27r) 3 P(k 1 )S^\J2 k -) * ( 22 ) 

a— 1 

3 

(C(k 1 )C(k 2 )C(k 3 )) = (2 7 r) 3 B(k 1 ,k 2 ,k 3 )^ 3 )(^k 0 ) , 

a= 1 

«(k 1 )C(k 2 )C(k 3 )C(k 4 )) 

4 

= (2^) 3 T(k 1 ,k 2 ,k 3 ,k 4 )^(^k 0 ). 


where <j>f denotes the value of the scalar field at the end 
moment of the phase tf when c = Co- The coefficient 
7 is an integration constant that can be determined by 
the background Friedmann equation. Without loss of 
generality, we keep 7 as a free coefficient and see its effect 
on the variation of the effective e-folding number. 

In our case, the VSL equivalent of the e-folding number 
is given by 

= H £)■ (27) 

and, correspondingly, the trajectories of the background 
dynamics can either be described by the phase space of 
or that of (AT, 7). In this regard, the method of 
calculating 6Af is similar to the analysis of j2l|. How¬ 
ever, the effects of the background evolution of the VSL 
cosmology upon curvature perturbations are analogous 
to those occurring in the study of the matter bounce 
cosmology [22l . [23|, where there is no quasi de Sitter ex¬ 
pansion. To be explicit, from now on we switch to the 
phase space of (</>, 7). Consequently, the VSL equivalent 
of the e-folding number takes the form: 

A7 = Ar«.,A) = ^r2/. (28) 

37 
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By expanding the scalar field 0, and determining the 
integration constant 7 via 0 —► 0 + S 0 and 7 —>■ 7 + S'), 
we can determine the curvature perturbations order by 
order up to third order: 


Ci=A/>0, (29) 

C 2 = ~^‘,<t><t>^4 > *h 007 + —Af^^S') , 

C3 = g'AC,00000 + gA/",00.y00 7 4" ^ A/" J 0'y'y007 —.A/]777*^7 ■ 

In the above expression, the subscript ^ = 9/90 denotes 
the derivative with respect to 0. After having obtained 
the expression A/ - in (l28l) , we can explicitly determine the 
coefficients such as A/", 0 . The above expression automat¬ 
ically includes the assumption that the field fluctuations 
are described before the primordial era by a highly Gaus¬ 
sian distribution. 

From ESI) we get 


7 = 7(0, 0) 


0 - 0 / 



(30) 


where W is the Lambert function. As a consequence, we 
can find 


h - = ,7 , H, (31) 

7 + 0-0/ 

where the contribution of <50 is secondary at superhorizon 
scales. Inserting the above expression into Eq. E9l) . we 
can derive 


^ 3(7 + 0 - 0 /)^’ 

^ 3(7 + 0 - 0 /) 1 2 ^ 

(27-0 + 0/) 3 

97(7 + 0 — 0 / ) 3 4 5 6 


(32) 

(33) 

(34) 


Recall that the nonlinearity parameters satisfy the re¬ 
lations in EU1) . Consequently, these parameters can be 
identified in the VSL cosmology as follows: 


/nl = 5, 

25(27-0 + 0/) 
Snl =-5- 


(35) 


From Eq. (1351) , we observe that the nonlinearity param¬ 
eter at second order /nl is a positive constant of order 
unity and thus is almost scale invariant. However, the 
nonlinearity parameter at third order <?nl is a function 
of the scalar field during the primordial phase and there¬ 
fore implies a strong scale dependence. 


IV. CONCLUSIONS AND DISCUSSION 


We have derived in our variable speed of light (VSL) 
model the primordial nonlinear fluctuations. We find 
that the parameter f^L — 5 which is compatible with 
the Planck2015 result: /nl = 2.7 ± 5.8 [l(|. The param¬ 
eter gpfL evolves during the primordial era and displays a 
running behavior. These result will allow the VSL model 
to be distinguished from inflationary models and other 
models such as the bouncing and string gas cosmology 
models. 
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